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Abstract 

A systematic investigation on muon anomalous magnetic moment and related lepton flavor- 
violating process such as /i — ► &y, r — > e^f and r — > is made at two loop level in the models 
with flavor-changing scalar interactions. The two loop diagrams with double scalar exchanges are 
studied and their contributions are found to be compatible with the ones from Barr-Zee diagram. 
By comparing with the latest data, the allowed ranges for the relevant Yukawa couplings Yj in 
lepton sector are obtained. The results show a hierarchical structure of Y^ e ,re ^ Yp T ~ Y m in the 
physical basis if Aa M is found to be > 50 x 10 -11 . It deviates from the widely used ansatz in which 
the off diagonal elements are proportional to the square root of the products of related fermion 
masses. An alternative Yukawa coupling matrix in the lepton sector is suggested to understand 
the current data. With such a reasonable Yukawa coupling ansatz, the decay rate of r — > [aj is 
found to be near the current experiment upper bound. 
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I. INTRODUCTION 



Recently, the Muon g — 2 Collaboration at BNL reported their improved result on the 
measurement of muon anomalous magnetic moment ( g — 2 ) [1] . Combining with the early 
measurements in CERN and BNL, the new average value of muon g — 2 is as follows 

a e * p = (116592030 ± 80) x 10~ u (1) 

This result confirmed the earlier measurement |2j with a much higher precision. With this 
new result the difference between experiment and the Standard Model (SM) prediction is 
enlarged again. The most recent analysis by different groups are given by 

( (303.3 ± 106.9) x 10- 11 [3] 
Aa M = a e * p - al M = I (297.0 ± 107.2) x 10- n (ex) [4] (2) 

{ (357.2 ± 106.4) x 10 _11 (m) [4] 

As the large muon g — 2 may imply the existence of new physics beyond the SM, in the 
recent years large amount of work has been done in checking the new physics contributions 
to muon g — 2 by using model dependent 0, 0] and independent approaches Q. 

In this work, we would like to focus on a general discussion on the models with lepton 
flavor- changing scalar interactions where the new physics contributions mainly arise from 
additional Yukawa couplings. Such models may be considered as the simple extension of the 
standard model (SM) with more than one Higgs doublet 0« {i > 1) but without imposing anv 
discrete symmetry. For example the extension of SM with two Higgs doublets (S2HDM)[8] 
motivated from the spontaneous CP violation [9!. Eo|. 

The general form of Yukawa interaction reads 

Cy = fcXSlfoa, (3) 

where Y°j (i, j = 1, 2, 3) is the Yukawa coupling matrices. The index a = 1, 2, ■ ■ • labels the 
Higgs doublets. The behavior of the Yukawa interactions depends on the texture of Yukawa 
coupling matrices. In general there are two kind of ansatz on Yukawa coupling matrices 
in mass eigenstates: 1. Yukawa coupling matrices of the scalar interactions are diagonal 
due to some discrete global symmetry [llj. 2. Yukawa coupling matrices contain non-zero 
off-diagonal elements which are naturally suppressed by the light quark masses |l2t , 13]. In 
the following sections (section II and IE) we discuss at two loop level the constraints on 
those Yukawa coupling matrix elements under the above two ansatz and will mainly focus 
on the latter one. One kind of two loop diagrams with double scalar exchanges are studied 
in detail and their contributions to muon g — 2 are found to be compatible with the one from 
Barr-Zee diagram. In section IV, combined constraints from muon g — 2 and several lepton 
flavor violating (LFV) processes are obtained. We note that unlike other experiments which 
often impose upper bounds of parameters in the new physics models, the current data on 
muon g — 2 may provide nontrivial lower bounds. It is founded that a small lower bound of 
Aa M > 50 x 10~ n will significantly modify the texture of Yukawa coupling matrix and make 
it deviate from the widely used ansatz in which the off diagonal elements are proportional 
to the square root of the products of related fermion masses. 



II. MUON g-2 FROM DIAGONAL YUKAWA COUPLINGS 

The ansatz of zero off-diagonal matrix elements is often used to avoid the flavor-changing 
neutral current (FCNC) at tree level which was originally suggested from the kaon physics, 
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such as K —y decay and K° — K mixing. Such a texture structure of the Yukawa 

couplings can be obtained by imposing some kind of discrete symmetries |l 11] . The minimal 
SUSY Standard Model (MSSM) and the two-Higgs-doublet model (2HDM) of type I and 
H can be cataloged into this type. In such models, the Yukawa interactions are flavor 
conserving and the couplings are proportional to the related fermion masses 

Y u = -f^-^i and V,, (I. (/ / ./) (4) 
Zrriw 

where g is the weak coupling constant and mw is the mass of W boson. £j is the rescaled 
coupling constant. In the minimal SUSY model and the 2HDM of type II, £j = tan/3(cot/3) 
for down (up) type fermions. The corresponding Feynman diagrams contributing to muon 
g — 2 at one loop level which is shown in Fig li b which have been recently discussed and 
compared with the current data in Ref. 




FIG. 1: One loop diagram contribution to muon g — 2 . The dashed curves represent the scalar or 
pseudo-scalar propagator, (a) Flavor conserving Yukawa interactions, (b) Flavor changing Yukawa 
interactions in which fx changes into r in the loop. 



As the muon lepton mass is small, i.e., <C m^, where is the mass of scalar (0 = h) 
or pseudo-scalar (0 = A) , the one loop contribution to muon g — 2 can be written as 0] 

87r 2 m2 \m*) 

where the sign "+ (-)" is for scalar {(f) = h) (pseudo-scalar <p = A) exchanges. It can be 
seen from the above equation that the one loop scalar contribution is not large enough to 
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explain the current data. Even for a large value of £ M = tan/? ~ 50, one still needs a very 
light mass of the scalar Mh ~ 5 GeV, which seems not favored by the LEP experiment. The 
situation will be even worse when both the scalar and pseudo-scalar are included as their 
contributions have opposite signs. 

The situation ma y b e quite different if one goes to two loop level. From the well known 
Barr-Zee mechanism |l8[] ( see. Fig|2])in which the scalar or pseudo-scalar couples to a heavy 
fermion loop. As the Yukawa couplings are no longer suppressed by the light fermion mass, 
the two loop contributions could be considerable. Taking the top quark loop as an example, 
the two loop Barr-Zee diagram contribution to muon (7 — 2 is given by 

h _ N c q 2 t m^m t ( m\ 
A % ~ 7J-2 m 2 F \y 4 ) Y«Y„ 

where N c = 3 and q t = 2/3 are the color number and the charge of top quark respectively. 
The integral function F(z) has the following formflil] 

— I L dx l ~^ x<yl r xS) In x ^~ x ^ for scalar 

2 JO x(l—x)—z z 

I Jq 1 dx x ^ x j_ z In x ^-~ x > for pseudo-scalar 

(7) 




FIG. 2: Two loop Barr-Zee diagram contribution to muon g — 2 



It is noticed that the contributions from Barr-Zee diagram through scalar and pseudo- 
scalar exchanges have also different sign, negative for scalar and positive for pseudo-scalar, 
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FIG. 3: Contribution to muon g — 2 from the two loop Barr-Zee diagrams. The three solid curves 
(from down to up) correspond to Y^(^) = 2 x 1(T 2 (48.3), 4 x 1(T 2 (97.6) and 1 x lO" 1 (273.9) 
respectively. The horizontal lines represent the la allowed rang from Ref.[3j 



which is just opposite to the one loop case. Thus there exists a cancellation between one 
and two loop diagram contributions. It was found in Refs.jl^, H3] that the pseudo-scalar 
exchanging Barr-Zee diagram can overwhelm its negative one loop contributions and results 
in a positive contribution to g— 2 . For a sufficient large value of the coupling £ M = tan (3 ~ 50, 
its contribution can reach the 2a experimental bound with < 70GeV. To avoid the 
cancellation between scalar and pseudo-scalar exchange, the mass of the scalar boson has 
to be pushed to be very heavy ( typically greater than 500 GeV). In FigJ3]the numerical 
calculation of Barr-Zee diagram contribution to muon ^7 — 2 is presented, which agrees with 
those results. 



III. MUON g-2 FROM OFF-DIAGONAL YUKAWA COUPLINGS 

When imposing the strict discrete symmetries to Yukawa interaction, the off-diagonal 
elements of Yukawa coupling matrix are all zeros. This is the simplest way to prevent the 
theory from tree level FCNC. However, to meet the constraints from the data on K° — K 
mixing and K — > the off-diagonal elements do not necessarily to be zero. An alterna- 

tive way is to impose some approximate symmetries such as global family symmetry |8( on 
the Lagrangian. This results in the second ansatz of the Yukawa matrices in which small 
off-diagonal matrix elements are allowed, which leads to an enhancement for many flavor 
changing processes. As the constraints from K° — K mixing are strong, the corresponding 
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off-diagonal matrix elements should be very small. However, up to now there is no such 
strong experimental constraints on the FCNC processes involving heavier flavors such as c 
and b quarks. The possibility of off-diagonal elements associated with the second and the 
third generation fermions are not excluded. One of the widely used ansatz of the Yukawa 
matrix basing on the hierarchical fermion mass spectrum m u ^ <C m c s <C m t ^ was proposed 
by Cheng and Sher [lil . ljj ]. In this ansatz, the off-diagonal matrix element has the following 



form: 



YiJ = & (8) 



w 



where are the rescaled Yukawa couplings which are roughly of the same order of mag- 
nitudes for all i,js. In this ansatz, the scalar or pseudo-scalar mediating d — s transition is 
strongly suppressed by small factor ^/mdm s / (2mw) — 4 x 1CT 4 , which easily satisfies the 
constraints from Am^,ex and T(K — > As the couplings grow larger for heavier 

fermions, the tree level FCNC processes may give considerable contributions in B° — B° 
mixing, fi + fi~ — > tc, \it and several rare B and r decay modes 21, 2^, IS H ■ 



Unlike the flavor-conserving one loop diagrams, the flavor-changing one loop dia- 
grams (see Fig^) with internal heavy fermions can give large contribution to muon 
g — 2 . The reason is that the loop integration yield an enhancement factor of ~ 
mjln(m?/m|)/(m At ln(m^/m^,)). For the internal r loop, it is a factor of 0(10). If one 
uses the scaled coupling ^ T and takes £ Mr ~ £ M as in the "Cheng-Sher" ansatz, the value of 
the enhancement factor can reach O(10 2 ). In the following discussion, for simplicity we only 
take one loop diagram with internal r loop into consideration as it is dominated over other 
fermion loops in the case that the Yukawa couplings are of the same order of magnitudes. 

The expression of one loop flavor changing diagram contribution to muon g — 2 is given 
by0 

1 m^m r ( m\ 3\ 2 



Aa„ = ±^-^1 In -4 - - Y' (9) 

M 8vr 2 mj V m l 2) ^ 

where the sign "+ (-)" is for scalar (<fi = h) (pseudo-scalar <p = A) exchanges. For detailed 
discussion on one loop flavor changing diagram, we refer to the Refs. |26L |27| ] . 

As the two loop contribution to muon g — 2 is more considerable via the Barr-Zee mech- 
anism in flavor-conserving case, it is nature to go further to consider the same diagram with 
flavor- changing couplings. However, in the case of muon g — 2 , as the initial and final states 
are all muons, it is easy to see that the Barr-Zee diagram with flavor- changing coupling can 
not contribute. It can only appear in flavor changing process such as \x — > erf. The non- 
trivial two loop diagrams which give non-negligible contribution to g — 2 are those diagrams 
(as shown in FigJHJ) which have two internal scalars with both of them coupling to a heavy 
fermion loop. 

It is known that large Yukawa couplings between scalar and heavy fermions can compen- 
sate the loop suppressing factor g 2 /167r 2 and make the Barr-Zee diagram to be sizable. The 
same mechanism also enhances the two-loop double scalar exchanging diagrams. Further 
more, in the flavor changing case, the \x lepton can go into heavier lepton r in the lower loop, 
this may provide an additional enhancement in loop integration. Taking the internal t-quark 
loop as an example, the ratio between the contribution to muon g — 2 from two- loop double 
scalar diagrams relative to the one from Barr-Zee type diagrams can be roughly estimated 
by the ratio between the couplings, which gives ~ ^^ T mim T /(4( (1 m^ sin 2 9w), where 6w 
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is Weinberg angle with the value sin 2 9w — 0.23. For the typical values of £t = 1 and 
£ MT ~ £ M = 30 the ratio is of order 1. Thus this kind of two- loop double scalar-exchanging 
diagram is compatible with the one of Barr-Zee type. In the large m t limit, the contribution 
to muon g — 2 from two loop double scalar (pseudo-scalar) exchanging diagram has the 
following form 

(10) 

The details of the two loop calculations can be found in the appendix at the end of this paper. 
Comparing with the one-loop flavor-changing diagram in the same way, one can see that 
the contribution from this diagram could be sizable . For a comparison, the contribution to 
muon g — 2 from two-loop double pseudo scalar-exchanging diagrams and Barr-Zee diagrams 
with pseudo-scalar are shown in Fig0] 



Aa^x 10") 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

Y 



FIG. 4: Comparison between two loop Barr-Zee pseudo-scalar and double pseudo-scalar exchanging 
diagram in contribution to muon g — 2 . The contribution to muon g — 2 is plotted as function 
of Y = Y^t = Yau- Three solid curves (from up to down) correspond to double scalar exchanging 
diagram contribution with scalar (pseudo-scalar) mass itla = 100, 150, 200GeV respectively. Three 
dashed curves indicate the ones from two loop Barr-Zee diagrams with pseudo-scalar exchange. 
The horizontal lines represent the ler allowed range from Ref.!3] 

To make the two kind of contributions comparable, we take = Y^ T = Y. It can 
be seen that the contribution from the former highly depends on the coupling Y and the 
scalar mass. In the range 0.05 < Y < 0.15, the contribution from double scalar-exchanging 
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diagram is much large than the one from Barr-Zee diagram when tua is about 100 ~ 150GeV. 
it decrease with itia increasing and becomes quite small when itla ~ 200 GeV. In Fig|3 the 
contribution to muon g — 2 from two-loop double scalar-exchanging diagrams is compared 
with the one from the corresponding flavor- changing one loop diagrams. Note that just 
like the case of Barr-Zee diagram, the two- loop double scalar- (pseudo-scalar) exchanging 
diagrams give negative (positive) contributions, which have the opposite signs as the one 
from one loop scalar (pseudo-scalar)- exchanging diagram. The reason is that a closed 
fermion loop always contributes a minus sign. It results in a strong cancellation between one 
and two loop diagram contributions in the case of flavor changing couplings with real Yukawa 
coupling constants. The allowed range of the scalar mass will be strongly constrained. 
Taking Y" MT = 0.08(^ T ~ 50) , Y tt = 0.67(& ~ 1) and Aa^ > 50 x 10~ n as an example, the 
mass of scalar lies in a narrow window of ~ 100 < < 200 GeV. 




FIG. 5: Comparison between one loop and two loop double scalar exchange diagrams in contribu- 
tion to muon g — 2 . The contribution to muon g — 2 is plotted as function of scalar mass. The two 
dashed curves represent the contribution at one loop with Y^ T (£ MT )=0.12(70.6) (up) and 0.08(47) 
(down) respectively. The two dotted curves correspond to the one from two loop double scalar 
diagram with the same couplings. (Note that their contribution are negative) The solid curves 
are the total contribution to g — 2 from the both diagrams. The horizontal lines represent the la 
allowed rang from Ref.0] 
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IV. LEPTON FLAVOR VIOLATION PROCESSES AND THE TEXTURE OF 
YUKAWA MATRIX 



The flavor changing Yukawa couplings will unavoidably lead to the enhancement of decay 
rates of lepton flavor violating processes, such as \i — > erf, r — > /i(e)7, \i — > e~e~e + and 
t — > e~e~e + (/i~/i~/x + ). The current experimental data especially the data of /J —>■ ey will 
impose strongest constraints on the related Yukawa couplings. From the current data the 
upper bound of the decay \i — > e-y is — > ej) < 3.6 x 10~ 30 GeV |28j|. It constrains the 
coupling Y er ( At ) to be extremely small. In the models with flavor changing scalar interactions, 
The leading contributions to \l — ► cy come from the one loop flavor changing diagram, 
the two loop double scalar exchanging diagram and the two loop flavor changing Barr-Zee 
diagrams. 

The effective vertex for one loop flavor changing scalar interaction reads 

rr = fa 4 - f) (ii) 

while the one for two loop loop double scalar exchanging is 



rr=^T?M-')W^ («) 

p ilit^m, \ mi 



In FigEl the decay rate from the sum of the first two diagrams are presented as function 
of the scalar mass. In the calculation we take the value of coupling Y tt = 0.67 (or £ t ~ 1). 
The value of Y^ T is taken to be 0.08(or £ MT ~ 50) which is the typical allowed value from 
the current data on g — 2 . It can be seen from the figure that the decay rate of \x — > ey 
constrain the value of Y Te to be no more than 10 -6 ~ 10 -5 for 100 < < 200GeV. 

Similarly, the value of coupling Y^ e is also constrained to be very small by the decay rate 
fi — > erf. The reason is that Y^ e is associated with the flavor changing Barr-Zee diagram in 
which muon goes into tau in the lower loop. If there is no accidental cancellation with other 
diagrams the upper bound of Y^ e can be obtained by assuming that the flavor changing 
Barr-Zee diagram is dominant. The decay rate of /x — > ey from this diagram alone can be 
obtained from Eq.fJUJ) and is given by 



T BZ (fx -> erf) = 8aml 



N c ql m^rrit 




Y a Y. 



(13) 



The numerical result is represented in Fig |7| which shows that the upper bound of Y^ e is also 
of the order 10" 6 ~ 10" 5 for 100 <m h < 200GeV. 

With the above constraints on the values of Yukawa couplings in the lepton sector, let us 
discuss the possible texture of Yukawa coupling matrix. In the SM with one Higgs doublet, 



it is well known that by assuming the Yukawa matrix to be of the Fritzsch form |30j, 131 
flavor basis, i.e. 



in 



v /m 1 m 2 
Y ~ ( y / m 1 m 2 v /m 2 m 3 | (14) 
v /m 2 m 3 m 3 



9 



r (n^ey) 




FIG. 6: Contribution to decay — > e^f from the sum of one loop and two loop double scalar 
diagrams. The three solid curves (from down to up) correspond to Y Te = 1 x 10~ 6 ,3 x 10~ 6 and 
1 x 10 -5 respectively. The coupling Y^ T is taken to be 0.08. The horizontal line indicates the 
experimental upper bound of \i — > e'y 



one can reproduce not only correct quark masses in mass eigenstates but also, in a good 
approximation, some of the mixing angles. In the models with multi-Higgs doublets, one 
can simply extend this Fritzsch parameterization to all the other Yukawa matrices including 
the leptons This results in the ansatz as in Eq.(jHJ) with all being of the same order 
of magnitude. 

It is not difficult to see that such an ansatz may be challenged by current experiment 
data in the lepton sector. This is because in order to explain the possible large muon g — 2 
, the off-diagonal elements connecting the second and third families should be enhanced, 
while to meet the constraints from /i — > ej, the ones connecting the first and second or the 
first and third families should be greatly suppressed. 



Taking the value of Aa M > 50 x 10 



-ii 



150GeV and ^> m/, as an example 



in the case of muon g — 2 , if the flavor- conserving Barr-Zee Diagram is playing the major 
role, the rescaled coupling £^ should be as large as 50 ( see FigOJ). If one assumes that the 
flavor-changing coupling is responsible for the large muon g — 2 , £ MT (Y" MT ) should be about 
10(0.02). On the other hand, due to the strong constraint from fi — > ej, for the flavor- 
changing contribution dominated case, £ Te has to be less than 0.08 when ^^ T {Y^ T ) is taken a 
typical value of 17.6(0.03). In the case of flavor changing Barr-Zee diagram dominant, the 
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r (n^ey) 




FIG. 7: Contribution to decay \x — > e7 from the two loop flavor changing Barr-Zee diagrams. 
The three solid curves (from down to up) correspond to Y^ e = 3 x 10 _6 ,7 x 10 -6 and 1 x 10 -5 
respectively. The horizontal line indicates the experimental upper bound of /i — > ej 



Yukawa coupling £ Me has to be less than 0.24. Thus one finds that 

^ ~ ^ T ~ O(10) 

U ~ ^ ^(lO- 1 ) (15) 

which clearly indicates that the rescaled couplings ^are not in the same order of magnitude. 
In the case of light pseudo-scalar mass rriA — 150 GeV and rrih ^> vtia the results are similar. 

From these considerations, it is suggested that the Yukawa matrices associated with the 
physical scalar bosons may take the following form in the mass eigenstate 

/ 0(1) 0(A n ) 0(A n ) \ 
Y ~l 2 I 0(X n ) 0(1) 0(1) (16) 
\0(A n ) 0(1) 0(1) / 

where / ~ 0.22 is roughly of the same order of the Wolfenstein parameter I, and n ~ 2 ~ 3. 
With such a parameterization, one is able to understand all the current experimental data 
concerning both muon g — 2 and lepton flavor-changing processes. 

If one takes the Yukawa matrix of the form in Eq.(|16|). the decay rate of r — > /ry could be 
predicted. In a good approximation, the decay rate can be obtained by replacing Y^ T Y Te into 
Y TT Y T/1 in Eqs. ljllj) and (|T2"j) Assuming r lepton dominance in the loop, the contributions to 
t — ► fi'y and shown in Fig|Hl The current upper bound on r — > fij is 3.5 x 10 -19 GeV[28j. 
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FIG. 8: Prediction of decay rate r — ► fxj from the sum of one loop and two loop double scalar 
diagrams. The three solid curves (from down to up) correspond to Y TT : 0.003, 0.01 and 0.03 
respectively. The coupling Y^ T is taken to be 0.08. The horizontal line indicates the experimental 
upper bound of r — > ^7 



It is found that, the predicted decay rate could reach the current experimental bound. A 
modest improvement in the precision of the present experiment for r — >• /ry may yield a first 
evidence of lepton family number non-conservation. 

In summary, we have studied the muon g — 2 and several lepton flavor violation processes 
in the models with flavor- changing scalar interactions. The two loop diagrams with double 
scalar exchanges have been investigated and their contribution to muon g — 2 is found to 
be compatible with the one from Barr-Zee diagram. The constraints on Yukawa coupling 
constants have been resulted from the current data of muon g — 2 and several lepton flavor 
violation processes. The results have shown a very strong constraints on the flavor-changing 
couplings associated with the first generation lepton. The early ansatz that the flavor 
changing couplings are proportional to the square root of the products of related fermion 
masses may not be suitable for the lepton sector if the Aa M is found to be > 50 x 10~ n . This 
indicates that both experimental and theoretical uncertainties need to be further reduced 
in order to explore the existence of new physics from muon g — 2 . It has been shown that 
an alternative simple parameterization given in Eq. (jl6|) is more attractive to understand 
the current experimental data. With such a parameterization, the decay rate of r —>■ fij is 
found to be close to the current experiment upper bound. 
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APPENDIX A: TWO LOOP DOUBLE SCALAR DIAGRAMS IN MUON g - 2 

AND n —> e7 

From the Yukawa interaction shown in Eq.(jnj). the qq<fi vertex has the following form in 
d dimension. 

ig^ 2 (Y 1 + Y 2l5 ) (Al) 

where \x is renormalization scale and e/2 = 2 — d/2 . The total amplitude can be written as 
the product of lower and upper parts as follows 

r„ = M ■ /„ (A2) 

The amplitude M for upper loop is given by 




FIG. 9: Two loop double scalar exchanging digram 



M 



-g 2 fi £ ■ 2N, 



c 



d d l 



(2 



7T 



[{A x + A 2l5 )(/l+ jq + m f )(B 1 + B 2l5 )(/l + m f )\ 



(I + q) 2 — m 2 . 1 2 — m 2 
(A3) 
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where Nc and m/ are the color number and mass of fermion /. For t quark / = t and 
Nc = 3. Ai 7 2 and -81,2 are the couplings of vertex A and B. 
The amplitude 1^ for the lower loop is given by 



J„ = -g 2 /2 e e(p 2 )(C 1 + C 2 i 5 )(h- A + mrhvUn- A + mr)(D 1 + D 2l5 )£(p 1 ) 
1 1 1 



{q — P2) 2 — m* (q-pi) 2 -m 2 (q 2 - mV) 2 



(A4) 



where C± :2 and D± i2 are the couplings for vertex C and D. 

After integrating over the lower loop and isolating the poles from Feynman integration, 
we obtain 



-8 • aWMA - A 2 B 2)mTl n, x r dx 2x(l _ f dy r-\ z{y + z){l _ y _ z) 



Jo 



(J^ - 2 lE + 21n47r - lnx(l -x) + ^)- h, div + 2 • / 1;Con ^ 

+ ^C ab R ■ ^ - 2 lE + 2 ln4vr - lnx(l - x)) • f 2 , div + 2 • / 2)< 

t{C x D x + C 2J D 2 + (C ± D 2 + C 2 D x ) lh )—^l 

2m, , 



(A5) 



with A' = (y + z)m 2 T + (l-y-z))m 2 , R = m}/[x(l - x)] and C ab = 2A 2 B 2 /(A 1 B 1 - A 2 B 2 ). 

In large m f limit, i.e. mj » » m 2 , the functions fi, div ( C on) and f 2 , div ( C on) have the 
following forms 



fl, d iv 
fl,con 



R 

ZV 2 ' 
R 
2ZV 2 



2,div 



A' 2 



1 -ln( 



A'R, 

1*' 



2, con 



2A' 2 



1 + ln 



A'i? 



(A6) 



After the renormalization in MS scheme for the upper loop, one finds 



r„ = 



-8 • 2N c g A m T m^A 1 B 1 - A 2 B 2 ) f 1 
(4tt) 4 J 

3 A' 4- R 1 A' 

+ (~ In x(l - x) + - + F{x) + In — ) + 2 • / 1)Con 



dx 2x(l — x) dy dz (y + z){l — y — z) 



A' 2 
1 1 



A' 2 2 



y ^1 y y 3A' 2R 

C ab R(- In x (l-x) + - + F{x) + In — )) + -C ab R(2 ■ f 2 , con - -) + ' 



2A /2 



I(p 2 )(C 1 D 1 + C 2 D 2 + (C7iD 2 + C 2 L> 1 )75))^^^(pi) 

2m,, 



(AT) 



with 



•F(x) = In 



m 2 — x(l — x)m 2 
V 2 



(A8) 
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In the limit of mj >> » m 2 , The above equation can be simplified as 

Ncg^m^m^A^ + A 2 B 2 ) / 5 m| 
M ~ 167r 4 m^ V 2 + n m 2 T 

{C X D X + C 2 D 2 )i(p 2 ) l -^e( Pl ) + (dD 2 + C 2 D x )t{ V2 ) ia ^ lh t{ Vl )} (A9) 

Therefore its contribution to muon g — 2 is as follows 

NcQ 4 m T m u m 2 f ( 5 rn\\ 
Aa„ = ~ ^ 167r I m r J (-J + In (^ + A 2 B 2 )(C lDl + C 2 D 2 ) (A10) 

In the real coupling case, for scalar exchange, one has 

g A 1 = gB l = Y ff , gd = gD 1 = V (All) 
and others are zero. Similarly for pseudoscalar exchange the couplings are 

g A 2 =gB 2 = iY ff , g C 2 = gD 2 = iY ilT (A12) 
Therefore, the two loop double scalar (pseudo-scalar) diagram's contribution to Aa M is 

2 / - m 2 



Ncm T m^m^ / 5 



A ^^^ri-2 +ln ^j r ^ < A13 > 

For the decay ji — > ey , the effective vetex is 



r( ^e 7) = N c g 4 m T m IM m 2 (A 1 B 1 + A 2 B 2 ) M ^ 

16vr 4 m4 V 2 ^ m 2 ) 



+ C^ 2 )l( P2 )^£( Pl ) + (C^ 2 + c5A)€(ft)^^€(pi) 

(A14) 

where C{ and C 2 are the Yukawa couplings for re0 vetex. The decay rate is then given by 

YU -> e 7 ) = — - — V lerff^VI 2 (A15) 
VP 16^^' M 1 v ; 
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